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Abstract: Optimal phase estimation of a phase-squeezed quantum state of light has 
been recently shown to beat the coherent-state limit. Here, the estimation is made robust to 
uncertainties in underlying parameters using a robust fixed -interval smoother 
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1. Introduction 

Precise estimation and tracking of a randomly varying optical phase is key to applications like communication [1], and 
metrology [2]. Optimal smoothed estimation of a widely varying phase under the influence of a stochastic Ornstein- 
Uhlenbeck (OU) process for a phase-squeezed beam noticeably surpasses the maximum accuracy that can be obtained 
for coherent light by 15 ±4% [3]. But, the estimation is very sensitive to fluctuations in the parameters underlying 
the phase noise or squeezing, and is desired to be made robust to such uncertainties. We have already demonstrated 
in Ref. [4], superior phase estimation as compared to the optimal case for a coherent state [5], using a robust fixed- 
interval smoother where the uncertain system satisfies a certain integral quadratic constraint (IQC) [6]. Here, the robust 
smoothing technique is applied to the case of squeezed state in comparison to the optimal case considered in Ref. [3]. 

2. Model of Adaptive Phase Estimation in Ref. [3] 

2.1. Process Model 

The phase 0(r) of the continuous optical phase-squeezed beam is modulated with an OU noise process, such that 

0(f) = -A</)(r) + y^v(f), (1) 

where A^' is the correlation time of 0(f), K determines the magnitude of the phase variation and v(f) is a zero-mean 
white Gaussian noise with unity amplitude. 

2.2. Measurement Model 

The phase-modulated beam is measured by homodyne detection using a local oscillator, the phase of which is adapted 
with the filtered estimate 0/(f ), thereby yielding a normalized homodyne output current /(f), 

I{t)dt ~ 2\a\[^{t) - (j)f{t)]dt + yfR^,dW{t), R,,, = aje-''" + (1 - aj)e-^'''", (2) 

where | a | is the amplitude of the input phase-squeezed beam, and dW [t ) is Wiener noise arising from squeezed vacuum 
fluctuations. The parameter R^g is determined by the degree of squeezing (r,„ > 0) and anti-squeezing (r^ > r„j) and 
by a'j = ([0(f) — 0/(f)]^)- We use the measurement appropriately scaled as our measurement model, 

0(f) := l/y^[/(f) + 2|a|0/(f)] -2|a|/y^0(f) + to(f), (3) 

where (o{t) = ^^-^^ is another zero-mean white Gaussian noise with unity amplitude. 

3. Robust Model 

3.1. Uncertain System 

We employ the robust fixed-interval smoothing technique from Ref. [6] as used in Ref. [4]. We introduce uncertainties 
8i and ^2, respectively, in the parameters X and 2\a\/ \/Rsq, such that Eq. (2.5) in Ref. [6] takes the form: 

Process: ^ = {-X + \/k/s.iK)^ + ^/KV , Measurement: 9 ^ \2\a\/ \ Rsq + ^iK] ^ + (o, (4) 
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: A2K(j) + 0) are the uncertainty inputs. 



where z = K(j) is the uncertainty output, and w = Ai K(j) + v and v : 
5.2. Robust vs. Optimal Smoothers for the Uncertain System 

We use the method from Ref. [4] to compute and compare the errors of the robust and optimal smoothers as a function 
of 5\ and &2, for various values of ^ and with = 1 .9 x lO"* rad/s, A = 5.9 x lO"* rad/s, | = 1 .00 x lO^s" ' , r„, = 0.36, 
and rp = 0.59 as used in Ref. [3]. Due to the implicit dependence of Rsq and (7^, we compute the error CJ,^ upon running 
several iterations until matches upto the 6''' decimal place to that of the prior iteration in each case. Fig. 1 shows 
the comparison for 20% and 40% uncertainties. Clearly, the robust smoother performs much better than the optimal 
smoother as 5\ and/or 82 approach — 1 for all levels of uncertainty. 
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Fig. 1. Comparison of Smoothed Errors for Uncertain System: (a) = 0.2, (b) /i = 0.4. 



4. Conclusion 

This work extends the treatment of robust smoothing to the case of squeezed state. For Rsq = 1, it boils down to the 
coherent state case of Ref. [4] when A2 = 0. These results may be demonstrated experimentally as part of further work. 
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